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We investigate the spin Hall effect in ballistic chaotic quantum dots with spin-orbit coupling. We 
show that a longitudinal charge current can generate a pure transverse spin current. While this 
transverse spin current is generically nonzero for a fixed sample, we show that when the spin-orbit 
coupling time is large compared to the mean dwell time inside the dot, it fluctuates universally from 
sample to sample or upon variation of the chemical potential with a vanishing average. For a fixed 
sample configuration, the transverse spin current has a finite typical value ~ e 2 V/h, proportional to 
the longitudinal bias V on the sample, and corresponding to about one excess open channel for one 
of the two spin species. Our analytical results are in agreement with numerical results in a diffusive 
system [W. Ren et al, Phys. Rev. Lett. 97, 066603 (2006)] and are further confirmed by numerical 
simulation in a chaotic cavity. 

PACS numbers: 72.25.Dc, 73.23.-b, 85.75.-d 



Introduction. The novel and rapidly expanding field 
of spintronics is interested in the creation, manipula- 
tion, and detection of polarized or pure spin currents 0] . 
The conventional methods of doing spintronics are to 
use magnetic fields and/or ferromagnets as parts of the 
creation-manipulation-detection cycle, and to use the 
Zeeman coupling and the ferromagnetic-exchange inter- 
actions to induce the spin dependency of transport. More 
recently, ways to generate spin accumulations and spin 
currents based on the coupling of spin and orbital de- 
grees of freedom have been explored. Among these pro- 
posals, much attention has been focused on the spin Hall 
effect (SHE), where pure spin currents are generated by 
applied electric currents on spin-orbit (SO) coupled sys- 
tems. Originally proposed by Dyakonov and Perel 0], the 
idea was resurrected by Hirsch 0| and extended to crys- 
tal SO field (the intrinsic SHE) by Sinova et al. 0] and 
Murakami et al. 0|. The current agreement is that the 
SHE vanishes for bulk, fc-linear SO coupling for diffusive 
two-dimensional electrons Jp, 0] ■ This result is however 
specific to these systems |9|, and the SHE does not van- 
ish for impurity-generated SO coupling, two-dimensional 
hole systems with either Rashba or Dresselhaus SO cou- 
pling, and for finite-sized electronic systems 0] • These 
predictions have been, to some extent, confirmed by ex- 
perimental observations of edge spin accumulations in 
electron and hole ^lj systems, and electrical detec- 
tion of spin currents via ferromagnetic leads |l2j . 

Most investigations of the SHE to date focused on dis- 
ordered conductors with spin-orbit interaction, where the 
disorder-averaged spin Hall conductivity was calculated 
using either the Kubo formalism or a diffusion equation 
approach 000IE0, 0, 0E1II1 Few numerical 
works alternatively used the scattering approach to trans- 
port 0| to calculate the average spin Hall conductance 
of explicitly finite-sized samples connected to external 



electrodes. These investigations were however restricted 
to tight-binding Hamiltonians with no or weak disorder 
in simple geometries 0, 0, . The data of Ref. 
in particular suggest that diffusive samples with large 
enough SO coupling exhibit universal fluctuations of the 
spin Hall conductance G s h with rms[G S H] ~ 0.18e/47r. 
These numerical investigations call for an analytical the- 
ory of the SHE in mesoscopic systems. It is the purpose 
of this article to provide such a theory. 

We analytically investigate the DC spin Hall effect in 
mesoscopic cavities with SO coupling. We calculate both 
the ensemble-average and the fluctuations of the trans- 
verse spin current generated by a longitudinal charge 
current. Our ap proach is based on random matrix the- 
ory (RMT) [2fJ, and is valid for ballistic chaotic and 
mesoscopic diffusive systems at low temperature, in the 
limit when the spin-orbit coupling time is much shorter 
than the mean dwell time of the electrons in the cav- 
ity- r so "C Td W oii [1^. We show that while the transverse 
spin current is generically nonzero for a typical sample, 
its sign and amplitude fluctuate universally, from sam- 
ple to sample or upon variation of the chemical potential 
with a vanishing average. We find that for a typical bal- 
listic chaotic quantum dot, the transverse spin current 
corresponds to slightly less than one excess open channel 
for one of the two spin species. These analytical results 
are confirmed by numerical simulations for a stroboscopic 
model of a ballistic chaotic cavity. 

Scattering approach. We consider a ballistic chaotic 
quantum dot coupled to four external electrodes via ideal 
point contacts, each with Ni open channels (i = 1, . . .4). 
The geometry is sketched in Fig. ^ Spin-orbit coupling 
exists only inside the dot, and the electrochemical po- 
tentials in the electrodes are spin-independent. A bias 
voltage V is applied between the longitudinal electrodes 
labeled 1 and 2. The voltages V3 and V4 are set such that 
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Figure 1: Ballistic quantum dot connected to four electrodes. 
The longitudinal bias V induces a charge current through 
terminals 1 and 2, while the voltages V-j,4 are adjusted such 
that no charge current flows through the transverse leads 3 
and 4. Spin-orbit coupling is active only in the gray region. 

no net charge current flows through the transverse elec- 
trodes 3 and 4. We will focus on the magnitude of the 
spin current through electrodes 3 and 4, in the limit when 
the openings to the electrodes are small enough, and the 
spin-orbit coupling strong enough that r so <C Tdweii- 
We write the spin-resolved current through the i-th 



electrode as 

The spin-dependent transmission coefficients are ob- 
tained by summing over electrode channels 

T if = EEiwi 2 , (2) 

i.e. £ m .cr;n,£r' is the transmission amplitude for an electron 
initially in a spin state a' in channel n of electrode j to a 
spin state a in channel m of electrode i. The transmission 
amplitudes t are the elements of the 2 Nt X 2Nt scattering 
matrix S, with Nt = Ylt=i N%- 

(z) 

We are interested in the transverse spin currents Z ( = 
ij — //, i = 3,4, under the two constraints that (i) charge 
current vanishes in the transverse leads, I J + l\ = 0, 
i = 3, 4 and (ii) the charge current is conserved, 1\ = 
—1-2 = I. From Eq. ljT]). transport through the system is 
then described by the following equation 
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where the transverse voltages (in units of V) read 

1 T (° 4 ) (T W _ off) + ( 2 jy 4 - j WW) _ t 
2 
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and we defined the dimensionless currents / = e 2 V J /h. 
We introduced generalized transmission probabilities 



Tr[(f 



ran J u b mn\ 



H = Q,x,y,z, (5) 



where arc Pauli matrices (a^ is the identity matrix) 
and one traces over the spin degree of freedom. 

Random matrix theory. We calculate the aver- 
age and fluctuations of the transverse spin currents j| , 
fi = x,y,z within the framework of RMT. Accordingly, 
we replace the scattering matrix S by a random unitary 
matrix, which, in our case of a system with time reversal 



symmetry (absence of magnetic field) and totally broken 
spin rotational symmetry (strong spin-orbit coupling), 
has to be taken from the circular symplectic ensemble 
(CSE) H3, E3, EE H3| We rewrite the generalized trans- 



(4a) mission probabilities T-^ as a trace over S 
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Here, m and n are channel indices, while a and (3 are 
spin indices. The trace is taken over both set of indices. 

Averages, variances, and covariances of the generalized 
transmission probabilities © over the CSE can be cal- 
culated using the method of Ref. j2f|. For the average 
transmission probabilities, we find 

(jM) = ( NiNi - -NA,) , (7) 

while variances and covariances are given by 
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where 5^ = 7 



Cm) _ /<j(t J -)\ 



y y > »j ' 
Because the transverse potentials V34 are spin- 
independent, they are not correlated with 7^ . Addi- 
tionally taking Eq. J7J into account, one concludes that 
the average transverse spin current vanishes (i = 3,4), 

= - jtf) E < T i M) >^> = °- ( 9 ) 

3=3,4 

However, for a given sample at a fixed chemical potential 
will in general be finite. Wc thus calculate var [J^j. 
We first note that (V 3 ,i)/V = (N^ N 2 )/2(N 1 + N 2 )., and 
that var [V3.4] vanishes to leading order in the inverse 
number of channels. One thus has 

var [J^\ = \ V varp#>] - \ covarpff^] (10) 
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+ 2cov^\7^}(V 3 )(V i ). 
From Eq. JSJl it follows that 



var [J, 
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Eqs. © and (|llf> are our main results. They show that, 
while the average transverse spin current vanishes, it ex- 
hibits universal sample-to-sample fluctuations. The ori- 
gin of this universality is the same as for charge trans- 
port [2(|, and relies on the fact expressed in Eq. © 
that to leading order, spin-dependent transmission cor- 
relators do not scale with the number of channels. The 
spin current carried by a single typical sample is given 
by rms[t7j ] X e 2 V/h, and is thus of order e 2 V/h in the 
limit of large number of channels. In other words, for a 
given sample, one spin species has of order one more open 
transport channel than the other one. For a fully sym- 
metric configuration, Ni = N, the spin current fluctuates 
universally for large N, with rmsfJf] ~ (e 2 V/h)/^/32. 
This translates into universal fluctuations of the trans- 
verse spin conductance with rmsfGW] = (e/47r\/32) ~ 
0.18(e/47r) in agreement with Ref. 

Numerical simulation. In the diffusive setup of 
Ref. |l9j the universal regime is not very large and thus 
it is difficult to unambiguously identify it. We therefore 



present numerical simulations in chaotic cavities to fur- 
ther illustrate our analytical predictions © and i|ll|) . 

We model the electronic dynamics inside a chaotic 
ballistic cavity by the spin kicked rotator [26L |27| . a 
one-dimensional quantum map which we represent by a 
2M x 2M Floquet (time-evolution) matrix |2S| 

7 ur = (UUXU^U) lv , M' = 0,1,...,M-1, (12a) 
Uu> = <We-^ +io)2/M <7o, (12b) 
U w = M- 1/2 e~ l2 * ll '/ M a Q , (12c) 

X W =6we -W/^)V(2*l/M)_ (12d) 

The matrix size 2M = 2L/\p 3> 1 is given by twice the 
ratio of the linear system size to the Fermi wavelength. 
The matrix X with 

V(p) = K cospco + K so (a x sin2p + a z ship). (13) 

corresponds to free SO coupled motion interrupted peri- 
odically by kicks described by the matrix II, correspond- 
ing to scattering off the boundaries of the quantum dot. 
In this form, the model is time-reversal symmetric, and 
the parameter l ensures that no additional symmetry 
exists in the system. The map is classically chaotic for 
kicking strength K > 7.5, and K so is related to the SO 
coupling time r so (in units of the stroboscopic period) 
through t so = 32tt 2 /K 2 M 2 [27| . From (T^, we con- 
struct the quasienergy-dependent scattering matrix |29j 



5(e) = P[e~ ie - 7(1 - P T P)}- 1 7P T , 
with P a 2Nt X 2M projection matrix 



P ka , k 'l3 



S a0 if/s' = ZW, 
otherwise. 



(14) 



(15) 



The l {k) (k = 1,2,..., 2N T , labels the modes) give the 
position in phase space of the attached leads. The mean 
dwell time Tdweii (in units of the stroboscopic period) is 
given by Tdweii = M/Nt- At large enough SO coupling, 
this model has been shown to exhibit the universality 
of the CSE. We refer the reader to Ref. [i?! for further 
details on the model. 

Averages were performed over 35 values of K in the 
range 41 < K < 48, 25 values of e uniformly distributed 
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Figure 2: Average and variance of the transverse spin current 
vs. the number of modes. Top panel: longitudinally sym- 
metric configuration with Nx = N2 = 2N 3 = 2N 4 = 2JV; 
bottom panel: longitudinally asymmetric configuration with 
N 2 = Ni = 2Ni = 2N 3 = 2N. In both cases the total num- 
ber of modes Nt = 6iV. The solid (dashed) lines give the 
analytical prediction JjjJ [1)11^] for the mean (variance) of the 
spin currents. Empty diamonds correspond to circles 
to var [J3 ] and triangles to var 



in < e < 2tt, and 10 different lead positions l^ k K We 
set the strength of K so such that t so = Td W cii/1250, and 
fixed values of M = 640 and l = 0.2. 

Our numerical results are presented in Fig. [21 Two 
cases were considered, the longitudinally symmetric 
(Ni = N2) and asymmetric (Ni 7^ N%) configurations. 
In both cases, the numerical data fully confirm our pre- 
dictions that the average spin current vanishes and that 
the variance of the transverse spin current is universal, 
i.e. it docs not depend on N for large enough value of N. 
In the asymmetric case N4 = 2N3, the variance of the 
spin current in lead 4 is twice as big as in lead 3, giving 
further confirmation to Eq. 

Conclusion. We have calculated the average and 
mesoscopic fluctuations of the transverse spin current 
generated by a charge current through a chaotic quan- 
tum dot with SO coupling. We find that, from sample 
to sample, the spin current fluctuates universally around 
zero average. In particular, for a fully symmetric configu- 
ration Ni = N, this translates into universal fluctuations 
of the spin conductance with rms[G S H] = (e/47rv32) ~ 
0.18(e/47r). This analytically establishes the universality 
observed numerically in Ref. 0] . 

We thank C.W.J. Bccnakker for valuable comments 
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